Algebraic and exponential decoherence in dissipative many-particle systems 
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The interplay between dissipation and many-body quantum correlations gives rise to a wealth of 
novel phenomena. Here we investigate spin-1/2 chains with uniform local couplings to a Markovian 
environment using the time-dependent density matrix renormalization group (tDMRG). For the 
open XXZ model, we discover that the coherence decay is algebraic instead of exponential. This is 
due to a vanishing gap in the spectrum of the corresponding Liouville superoperator - irrespective 
of the Hamiltonian gap - and can be explained on the basis of a perturbative treatment. On the 
contrary, decoherence in the open transverse-field Ising model is found to be always exponential. In 
this case, the interactions can both facilitate and impede the environment-induced decoherence. 
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I. INTRODUCTION 

If quantum systems are well-isolated from their envi- 
ronment, their time evolution is unitary as described by 
the Schrodinger equation. Ultimately though, all quan- 
tum systems are inevitably coupled to their surround- 
ings which may fundamentally alter the static and non- 
equilibrium properties [1-3]. These open quantum sys- 
tems provide novel perspectives for quantum many-body 
physics where the interplay between the many-particle 
correlations and dissipation gives rise to a wealth of novel 
phenomena, challenging both theory [4-18] and experi- 
ment [19-21]. Among the features of open quantum sys- 
tems, the concept of decoherence, the dynamical destruc- 
tion of quantum coherence due to interaction with an en- 
vironment, is of paramount importance as it poses fun- 
damental limitations on quantum information processing 
as in quantum cryptography, communication, and com- 
putation. On the other hand, one can try to exploit 
the effect and engineer the dissipative processes to yield 
desired dark states which introduces a new concept for 
quantum state preparation [22-25]. 

Decoherence has been studied extensively for small or 
non-interacting quantum systems. For interacting dissi- 
pative many-particle systems however, the effect of the 
quantum many-body correlations on the decoherence is 
far from being fully understood. For a given system- 
environment coupling, it is, for example, not clear from 
the outset whether interactions in the system will fa- 
cilitate or impede decoherence. Recently, this prob- 
lem has been addressed for a two-site Bosc-Hubbard 
model, which can be mapped to a classical diffusion pro- 
cess exhibiting non-Brownian behavior and an imped- 
ing of the decoherence [11]. In this work, we use the 
time-dependent density matrix renormalization group 
(tDMRG) [26-28] to investigate the decoherence dynam- 
ics in genuine open quantum many-body systems and 
explain the results on the basis of an additional perturba- 
tive analysis. We focus on the problem of how the corre- 
lated many-body effects fundamentally change the deco- 
herence dynamics and give rise to novel phenomena that 



are absent for the corresponding small or non-interacting 
systems. 

Considering a spin-1/2 chain uniformly coupled to a 
Markovian environment, the time evolution of the system 
state p is governed by a Lindblad master equation [29] 
(we set h=l) 



d t p = Cp = lip + Vp 

= -i[H, p] + 7 E ( L ^ L l - \i L \ L ^ P} 



(1) 



Here, C denotes the Liouville superoperator containing 
two parts: T-Lp — —i[H,p] generates the evolution due to 
the system Hamiltonian H while the dissipative process is 
described by Vp. L t are the Lindblad operators describ- 
ing the coupling of the system to the environment and 
7 quantifies the strength of this coupling. Throughout 
the paper we choose Li = Sf . This type of coupling was 
first introduced in the study of dissipative two-state sys- 
tems [30] and is widely used to describe the environment- 
induced decoherence. A short discussion of experimen- 
tally and technologically relevant realizations is given in 
section IV. Consider a single qubit (spin-1/2) in the ab- 
sence of external fields (H — 0) that is coupled to the en- 
vironment via L\ = Sf. The master equation (1) predicts 
that the off-diagonal element of the qubit density matrix 
p decays exponentially with time, pf,i{t) ~ e~ 7 *' 4 , which 
means that superpositions of quantum states (quantum 
coherence) are rapidly destroyed. Below, we investigate 
the long-time dynamics of one or two neighboring spins 
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Figure 1: A quantum spin chain uniformly coupled to the 
environment via the z-components of the spins. 



in the center of an open spin- 1/2 chain. The interplay 
between interaction and dissipation qualitatively alters 
the decoherence behavior. In particular, we find (i) that, 
for the Heisenberg XXZ model, the coherence decay be- 
comes algebraic instead of exponential and (ii) that, for 
the transverse-field Ising model, the coherence decay re- 
mains exponential but that many-body interactions can 
both facilitate and impede the decoherence. We explain 
both features on the basis of a perturbative treatment for 
the Liouville superoperator. 

Before addressing the aforementioned specific models, 
some general remarks are appropriate. First, as long 
as L\ — Li Vj, it is easy to check that the maximally 
mixed state po °c 1 is always a steady state (d t p = 0) 
of Eq. (1). For the models addressed in this paper, p 
or restrictions of it to certain symmetry sectors are the 
unique steady states. Though all the initial states will 
eventually converge to such a steady state, the approach 
towards it is typically highly nontrivial and depends on 
the quantum many-body Hamiltonian. The dynamics 
(1) is governed by the non-Hcrmitian superoperator C. 
We denote its eigenvalues by A Q . They have non-positive 
real-parts, Re A Q < 0, and the steady state has the eigen- 
value Ao = 0. We call 
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min Re (— A Q ) 

a>0 ' 



(2) 



the spectral gap of the Liouville operator. It can be used 
to characterize the long-time behavior of the system. If 
the gap is finite, the distance of the time-evolved state to 
the steady state will decrease exponentially with time, 
and A sets the corresponding inverse relaxation time. 
However, as we will see below, the many-body interac- 
tions in the system may qualitatively alter the dynamics 
by closing the gap of C in the thermodynamic limit which 
gives rise to a novel algebraic long-time decoherence be- 
havior. 



II. ALGEBRAIC DECOHERENCE IN THE 
OPEN XXZ MODEL 



First, let us consider the one-dimensional (ID) 
1/2 XXZ model 
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(3) 



uniformly coupled to the environment via the z- 
components of the spins, L{ = Sf. We study the time 
evolution of the system density matrix p(t) based on the 
master equation (1) with the initial state p(0) = \^ ) (^ \ 
being the Neel state |\?o) = |"fl ••• ||). In the absence 
of dissipation (7 = 0) , the time evolution for this setup 
has for example been studied in the context of quantum 
quenches [31-33], where the long-time behavior highly 
depends on J z /J xy - In this model, the total magneti- 
zation J2i S* is conserved. As a consequence, the off- 
diagonal element pi , = (Sf) of the single-site density 
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Figure 2: (a) Power-law decay of the off-diagonal density ma- 
trix element C = {SfS~ +1 } in the dissipative XXZ chain (3) 
of length N — 96 for different J xy , J z and fixed bath coupling 
7 = 1, evaluated in the center of the chain, (b) Finite-size 
scaling of the gap (2) of the Liouville superoperator C for the 
open XXZ model, obtained by exact diagonalization (ED). 



matrices are strictly zero for all times and can not be 
used to monitor the decoherence. Instead, we choose the 
off-diagonal term C = {Sf S^ +1 ) of the two-site density 
matrix to quantify the decoherence, where sites i and i + 1 
are located in the center of the chain. For the simplest 
case of a two-site system (N = 2), it is easy to show that, 
in the subspace of zero magnetization, the model can be 
mapped to the decoherence problem of a single spin sub- 
ject to a transverse field and that the off-diagonal element 
of the reduced density matrix decays exponentially (see 
appendix B). 

In order to study the effects of the many-body cor- 
relations on the decoherence, we employ tDMRG [26- 
28]. As shown in Ref. [34], the propagator exp(£t) can 
be approximated by a circuit of two-site gates with an 
accuracy that is well-controlled in terms of the opera- 
tionally relevant (1 — > l)-norm. Here, we specifically em- 
ploy a fourth-order Trotter-Suzuki decomposition with 
a time-step of size At = 0.125. Starting from certain 
product states p(0), the time-evolved states are obtained 
by applying the local Trotter gates and approximat- 
ing pit) by matrix product states (MPS) [35-37]. The 
essence of the DMRG procedure is to express p(t) in ev- 
ery step of the simulation in a reduced Hilbert-Schmidt 
orthonormal operator basis {Of ®Of} for a spatial split- 
ting of the system into a left and a right part so that 
p{t) = J2i v iOf ® Of. Such a representation can always 
be obtained by singular value decomposition. The ap- 
proximation consists in discarding all components i with 
weights vf/J2j u j below a certain threshold e (between 



10~ 10 and 10~ 12 in this work). One has to ensure conver- 
gence of the numerical results with respect to the trunca- 
tion threshold e and the system size N in order to capture 
the physics of the thermodynamic limit. This is carried 
out in appendix C. 

The evolution of the coherence C(t) in open XXZ 
chains is shown in Fig. 2. In the long-time limit, we 
find the coherence to decay algebraically, instead of ex- 
ponentially, according to the power-law 



C{t) ex t' 



ith 



rj: 



1.58. 



(4) 



The exponent 77 is in the studied parameter regime in- 
dependent of the system parameters J xy , J z , and 7. In 
general, an algebraic decay implies the absence of a char- 
acteristic time scale in the long-time dynamics. It results 
from the vanishing of the gap A of the Liouvillian C in 
the thermodynamic limit. That this is indeed the case 
can be verified numerically by exact diagonalization as 
shown in Fig. 2b. 

To get a better understanding for the vanishing of 
the gap A, let us perform a second order perturba- 
tive analysis to derive an effective superexchange Li- 
ouvillian £ c ff for the limit 7 3> | J z \ , \ J xy | of strong 
dissipation. The Liouvillian can be split into an un- 
perturbed part Cqp := — i[H z ,p] + T>p, where H z = 

:= -i[H xy ,p] 

The steady 

<r)(cr|, where 



^Ei^f^+i' an d th e perturbation dp 



with H r „ — 



states of £0 (eigenvalue A = 0) are p° = 
\tr) = |(7i . . .ctat) are the {5f}-eigenstates spanning the 
Hilbert space of the spin configurations with zero total 
magnetization. The effect of a small coupling J xy is to 
lift the degeneracy in the steady state manifold through 
a superexchange process which leads us to an effective Li- 
ouvillian £ c ff, constrained to the subspace EI spanned by 
the operators p$ . H is to be understood as a subspace of 
the vector space B(J4?) of linear operators on the Hilbert 
space Jff, i.e., t : B(Jf?) -> B(Jt?) and £ off : H -> H. 
One obtains 



£ e ff = VCl 



1 



^o — £0 



-CxP. 



(5) 



V is the projector onto the subspace EL The interme- 
diate states in the perturbation theory are of the form 
Ar' - W)W'\, where \*>) = (S+Sr +1 + <?r<?+ J^) for 

some bond (i,i + 1). Their >C -eigenvalues, needed to 
evaluate the denominator in Eq. (5), are —7 or — 7±iJ z , 
depending on <x. However, the term ±iJ z can be ig- 
nored as it represents an irrelevant contribution of order 
I/7 2 to C e g. The full calculation given in appendix A 1 
shows that the matrix elements of the effective Liouvillian 
are identical with those of the ferromagnetic Heisenberg 
model 
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in the sense that 
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(7) 



As a consequence, at the level of the second order per- 
turbation theory, the gap (2) of the effective Liouvillian 
£ e fi is that of the Heisenberg ferromagnet H e g. Its gap 
vanishes as 1/A^ 2 due to the quadratic spin-wave disper- 
sion around zero momentum and the 2n/N spacing of the 
quasi-momenta. This explains the quadratic behavior of 
the gaps A for the full model in Fig. 2b. 
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Figure 3: Time evolution of |pi-j,(t)| for the dissipative 
transverse-field Ising model (8) with different J z , fixed bath 
coupling 7 = 1, field h x — 1, and system size iV = 64. 



III. DECOHERENCE IN THE OPEN 
TRANSVERSE ISING MODEL 

A second paradigmatic example is the dissipative 
transverse-field Ising model 

i 

= J2(4J*SlS? +1 -2h x S?) (8) 

i 

with the interaction strength J z , the transverse magnetic 
field h x , and the Pauli matrices af . To study the inter- 
play of interaction and dissipation, we set for simplicity 
h x = 1 and vary J z and the bath coupling 7. As Lindblad 
operators, we choose again L. t = Sf and study the time 
evolution of the system density matrix based on Eq. (1) 
starting from a fully polarized state, i.e., p(0) = l^gK^ol 
with |*o) = j TT * * * TT) - Alternative initial states have 
also been checked. However, as we are foremost inter- 
ested in the long-time behavior, the choice of the initial 
state is of minor importance. Let us first consider the 
nonintcracting case with J z = which reduces to the de- 
coherence problem of a single spin subject to an external 
field. In this case, the off-diagonal element pi , = (S^) 
of the single-site reduced density matrix decays exponen- 
tially in the long-time limit as p\ At) ~ e~ A °*, where the 
decay rate is A = 7/4 (as long as 7 < 8\h x \; see ap- 
pendix B). For the interacting many-body system, one 
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Figure 4: (a) The renormalized decoherence rate A/7 [Eq. (9)] 
for the dissipative transverse Ising model (8) as a function of 
J z for different bath couplings 7. (b) Time evolution of the 
magnetization (S*) for fixed 7 = 1 and different J z , where site 
i is in the middle of the chain, (c) Dynamics of (Si) for fixed 
J z = 0.25 and different 7. In all cases the system contained 
N = 64 sites. 



the dynamics of the magnetization (S z ) for different J z 
and 7, as shown in Figures 4b and 4c. For small J z and 7, 
(Sf) is well described by an exponentially decaying os- 
cillation. Generally, for open many-body systems, the 
long-range quantum correlations are usually destroyed 
during the long-time dissipative dynamics. As a conse- 
quence, the quantum entanglement between a single spin 
and the rest of the system is weak. This allows us to ex- 
plain the above observations in a mean-field framework 
corresponding to the decoupling of the interaction term 
SfS z +1 ~ S Z (S Z ). On a qualitative level, the decoher- 
ence in the long-time limit can be understood as that of 
a single spin in a constant transverse field and a longitu- 
dinal field 2(S z (t)} due to its nearest neighbors. Figures 
4b and 4c show that, for small J z and 7, the longitudinal 
field (S z (t)) is quickly oscillating - hence, playing a role 
similar to that of noise and thus accelerating the deco- 
herence. Once, the oscillations of (S z (t)) vanish (large 
Jz or 7), the decoherence is suppressed. 

The second key observation, that strong interaction 
impedes the decoherence, can again be understood on 
the basis of a perturbative analysis, here in the limit of 
a weak magnetic field, 7> \h x \. The field terms ex h x of 
the Liouvillian are considered as a perturbation so that 
C = d + d with 



can use \pl , (t)\, with site i in the middle of the chain, to 
monitor the coherence decay. In contrast to the situation 
for the open XXZ model, we find that in the long-time 
limit the coherence always decreases exponentially 



\P\,M 



-At 



(9) 



where the decoherence rate (inverse relaxation time) A 
is determined by the interplay of interaction and dissipa- 
tion. 

The evolution of \pi ,(t) | for different J z and fixed bath 
coupling 7 = 1 is shown in Fig. 3. For small J Zl the de- 
coherence dynamics is well described by oscillations of 
exponentially decaying amplitude. For large J z , \pi ,(t)\ 
decays exponentially without oscillations. Similar behav- 
ior was found in the relaxation dynamics of the quenched 
ID XXZ model without dissipation [31-33]. Now we turn 
to the question of whether the interaction facilitates or 
impedes the decoherence. In other words, the question is 
whether the decoherence rate A is below or above that 
of the noninteracting case J z — with Ao = 7/4. As 
shown in Fig. 4a, the answer to this question depends 
in an intricate manner on the values of 7 and J z . In the 
presence of weak dissipation (small 7), the dependence of 
A on J z is non-monotonic. The interaction facilitates the 
environment-induced decoherence (A > Ao) for small J 2 , 
while it impedes the decoherence (A < A ) for large J z . 
For sufficiently strong dissipation, the interaction always 
suppresses the decoherence. 

In order to obtain a better understanding of the fact 
that interaction and dissipation cooperate to enhance the 
decoherence in the case of small J z and 7, we calculate 



C p = -i[H z ,p]+T>p and dp = -i[H x , p], (10) 

where H z = 4J, E< Sf 5f +1 and H x = -2/i^Sf. In 

the second order perturbation theory, the eigenoperators 
of d arc similar to those in the treatment of the open 
XXZ model (now, the dynamics is not constrained to 
sectors of constant magnetization), but the intermediate 

states are different. Their £ ~ e ig enva l ues are ~~ 7/2 and 
— -f/2±i4J z . The effective Liouville superoperator (5) is 
again of the form (7) and the effective Hamiltonian reads 



HqS = 



Z^ 



- (a - a^Sf^S;, 



i+1 



2 ** 



(11) 
in this case, where a — 16h x /j and a' — Ah x r y/('-{ 2 /A + 
(4J Z ) 2 ). Based on Eq. (11), one can show that the gap 
(2) of the effective Liouvillian has for small J z a value 
< (a + a')/4. For sufficiently large J z , the gap is given 
by a' . The corresponding eigenstate of H g is the spin- 
wave-like state J2 3 -I t x ••• f c i]Vj+it x ■■■ t x )> where 
<Sf|tf) — |ltf)- A detailed derivation is given in ap- 
pendix A 2. So the gap decays as 1/Jf, i.e., strong in- 
teraction impedes decoherence as we have found in the 
quasi-exact numerical analysis. 



IV. EXPERIMENTAL REALIZATIONS AND 
APPLICATIONS 

Besides being of fundamental theoretical interest, the 
two dissipative models addressed in this paper are of 
broad experimental and technological relevance. We 



shortly mention a few examples. A uniform coupling 
to the environment via Lindblad operators Li — S* oc- 
curs for example naturally in quantum computer archi- 
tectures [38, 39] based on superconducting flux qubits 
(rf-SQUIDs) through fluctuations of the external mag- 
netic flux [40-42] . Inductive coupling of flux qubits yields 
the Ising-type interaction SfSj [43-45]. In ultracold 
atom systems [46] where both interaction and dissipa- 
tion can be controlled, the corresponding Lindblad oper- 
ators Li — m describe laser fluctuations and incoherent 
scattering of the laser light [47-49] . With quantum dot 
spin-qubits [50, 51] one can implement both, the trans- 
verse Ising model [52] and the Heisenberg model [50], 
where Li = Sf describes the effect of variations in the 
longitudinal nuclear magnetic field. 



V. CONCLUSION 

In summary, we have studied the long-time dynam- 
ics of open quantum spin systems, discovering that the 
quantum many-body effects can significantly change the 
nature of the environment-induced decoherence by either 
altering the exponential coherence decay to being alge- 
braic or, alternatively, by increasing or decreasing the 
decay rate. Besides illustrative quasi-exact tDMRG sim- 
ulations, we have explained those effects by a pertur- 
bative analysis. Another interesting direction for further 
investigations are driven-dissipative quantum many-body 
systems, where external force drives the system far from 
equilibrium and the interplay between driving, dissipa- 
tion, and interaction may give rise to further novel non- 
equilibrium phenomena. 
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Appendix A: Perturbative treatment 

1. Open spin-1/2 XXZ model with Li = St 

Let us again consider the spin-1/2 XXZ model (3) cou- 
pled to a Markovian environment according to the Lind- 
blad master equation (1) with uniform Lindblad opera- 
tors Li = Sf and bath coupling strength 7. Based on 
a perturbative treatment, one can derive an effective Li- 
ouville superoperator £ e g- for the limit of strong bath 
coupling 7 ^> \J X y\, \Jz\- To this purpose, the Liouvillian 
is split into an unperturbed part Cq and a perturbation 



C\ such that 



t = to + t u (Al) 



C\p = -i[H xy ,p] 



(A3) 



where H z = J z £, SfSf +1 and H xy = % E t (S+Sr +1 + 
S, t S^ +1 ) as in Eq. (3). The Liouvillian conserves the 
total magnetization such that [^ Sf, p(t)] = V4. As 
in the main part of the paper, we choose the symmetry 
sector of zero magnetization for the following. 

The steady states (eigenvalue A = 0) of the unper- 
turbed part £0 are of the form p% := \er)(er\ with {S 1 ?}- 
eigenstates \<r) — \(J\...(7n) obeying Sf\cr) = ^f|er). 
Within the second order perturbation theory, the result- 
ing effective Liouvillian is given by Eq. (5) with V pro- 
jecting onto the space H, spanned by the unperturbed 
steady states p$. Starting from such a state, the inter- 
mediate states (after application of £1) are of the form 
A"' := |<r)<<r'| such that \*>) = (S+Sr ±1 + S~ S+ +1 )\a) 
for some bond (i, i + 1). We need their Co eigenvalues to 
evaluate the denominator in Eq. (5). Computing CoA a ' a ' , 
one finds that the eigenvalues arc —7 or —7 ± i J z de- 
pending on the values of c,_i, <7j (ctj+i = — Cj), and 
CTi+2- However, as 7 3> \J xy \, \Jz\ and one only needs the 
effective Liouvillian up to first order in 1/7, the factor 
l/(Ao — Cq) in Eq. (5) can for all cases be approximated 
by I/7. 



1 1 



7 ± Mz 7 



+ 0(T 2 ) 



With the terms h % := ^f- (Sf S~ +1 + S i Sf +1 ) in H xy , the 
effective Liouvillian attains the form 

-VUUft = --V[H xy , [H xy ,pZ]} 

7 7 

= -" E ( h iPo + Potf - IhipZhi) . (A4) 



Using that hf — -^ (| — SiSf +1 ), one can read off the 
effective Hamiltonian H c r that describes the action of the 
effective Liouvillian C c g on the space H according to the 
definition (7). It turns out to be the isotropic Heisenberg 
ferromagnet with coupling constant J xy /j, 



Hc.tr — — 



7 2 

u X y ST~^ 



\"i "i+1 + "» "i+l) + "» "»+l 



1" 

4 

(A5) 
The spectra of — C e g and H c fr are identical. The gap 
is decisive for the long-time dynamics and decoherence. 
To confirm and illustrate the perturbative analysis, we 
have done an exact diagonalization of the full model C 
and the effective model C c rr. As shown in Table I, with 
increasing 7, the spectrum of C e g converges indeed to- 
wards the corresponding eigenvalues of C. The exact gap 



of £ e ff (finite only for finite system sizes N) is available, 
because spin waves with the well-known cosine dispersion 
relation are exact eigcnstatcs of the isotropic Heisenberg 
fcrromagnct (A5). 



n 


-t 


—C«B 

7= 100 


— £eff 

7 = 10 














1 


0.292893 


0.292877 


0.291271 


2 


0.633975 


0.633925 


0.629092 


3 


1.000000 


0.999994 


0.999369 


4 


1.707107 


1.707042 


1.700617 


5 


2.366025 


2.365987 


2.362134 



Table I: Spectrum of the effective model £ e ft [Eq. (A5)] com- 
pared to the corresponding eigenvalues of the full Liouvillian 
[Eq. (Al)] for a XXZ chain of JV = 4 sites with open boundary 
conditions, J xy = 0.5, and J z — 1. The eigenvalues are given 
in units of the effective coupling J^y/l- 



2. Open transverse Ising model with Li = S* 

This appendix presents, in full detail, the perturbation 
theory for the spin-1/2 Ising model in a transverse field 
(8) coupled to a Markovian environment according to 
the Lindblad master equation (1) with Li = Sf and bath 
coupling strength 7. Using a second order perturbation 
theory, one can derive an effective Liouville superopera- 
tor £ e g for the limit of a weak magnetic field 7 ^> \h x \. 
To this purpose, the Liouvillian is split into an unper- 
turbed part £q and a perturbation C\ such that 

£ = £o+£i, (A6) 

C aP = -i[H z ,p]+^(SfpSf - P/*>> ( A? ) 



£ip = -i[H x ,p], 



(A8) 



where H z = 4J Z £. SfSf +1 and H x = -2h x £. Sf as in 

Eq. (8). 

As in the perturbative analysis of the XXZ model, the 
steady states (eigenvalue Ao = 0) of the unperturbed part 
£ are of the form p% := |<x)(<x| with {Sf }-eigenstates 
\a) obeying Sf\cr) — ^\cr). The effective Liouvillian 

is given by Eq. (5) with V projecting onto the space H 
that is spanned by the unperturbed steady states p$ . 
The differences to the calculation for the XXZ model in 
appendix A 1 are that, here, the total magnetization is 
not conserved and that the factor l/(Ao — £0) in Eq. (5) 
is now nontrivial. Starting from one of the unperturbed 
steady states, the intermediate states (after application 
of £1) are of the form A CTCT ' := |<x)(er'| such that \<r') = 
2Sf \tr) for some site i. Computing £oA CTCT , one finds 
that the £0 eigenvalues of the intermediate states are 
-7/2 - %{E a - E a ,). Here, E„ := (<r\H x \a) and E„ - 



Eu-t = or ±4</ z depending on the values of <7j_i, <7j, and 
Uj+1. With the definition of the state \tr,i) := 2Sf \a) 
and the corresponding eigenvalue of H z labeled by E^, 
the effective Liouvillian (5) takes the form 



£ cH \ ( t}(<t\=V£ 1 



1 



i 



AkXH 

h x \<T,i)(cr\ 
1 /2 + i(E] T -E (T ) 

\<t){(t\- \(T,i)((T,i\ 
( 7 /2) 2 + (E% - E a ) 



h.c. 



(A9) 



Employing |i^. — _E CT | = 2| J z {cfi-\ + o»+i)|, one can read 
off the effective Hamiltonian H s that describes the ac- 
tion of the effective Liouvillian C e s on the space EI ac- 
cording to the definition (7) and finds 



H e e = 



*' E ( A + S i-l S i+l ) \9~ S i 



E 



-{a-a')SUS. 



i+l 



o~ 



with 



7(2/i, 



(7/2) 



and a 



7(2/**) 



( 7 /2) 2 + (u z y 



(A10) 



(AH) 



We are particularly interested in the gap (2) of £ c h 
which is given by that of H c g and determines the deco- 
herence. Let |ff ) and \\.f) denote the normalized eigen- 
states of Sf with eigenvalues ±1/2. In this notation, the 
fully polarized state 



iVo) = irr---r) 



(A12) 



is the eigenstate of H e g with eigenvalue Ao = 0. As H c g is 
positive-semidefinite (H e g y 0), \tf>o) is its ground state. 
Note that the counterpart of \ipo) in H is the maximally 
mixed state 1 = J2& \ cr ) (""I which is also the exact steady 
state of £. 

In order to determine or bound the gap, the first intu- 
ition may be to flip a single spin, i.e., to consider states 



\4>\)-.= 1 rv ■■■vi*v 

They provide the upper bound 



r>. 



^I^effl^l) 



a + a 
~ 4 



(A13) 



(A14) 



to the gap. The same bound pertains for any superpo- 
sition of the states \4>\). One can also construct exact 
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Figure 5: Convergence analysis for the open XXZ model (Al) at the example of the expectation value (SfS~ +1 ) with different 
J xy and J z , and fixed bath coupling 7 = 1, where site i is in the middle of the chain. Left: variation of the DMRG truncation 
threshold e for the system size N = 96. Right: variation of the system size for the truncation threshold e = 10 -12 . 
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Figure 6: Convergence analysis for the open transverse-field fsing model (A6) at the example of the magnetization {Sf) with 
different J z , fixed bath coupling 7 = 1, and field h x — 1, where site i is in the middle of the chain. Left: variation of the DMRG 
truncation threshold e for the system size N = 64. Right: variation of the system size for the truncation threshold e = 1CP 10 . 



spin- wave- type eigenstates \ijj^ ) where a pair of spin de- 
fects is bound and travels through the lattice. 

i^) : = iff . . . ri*i*+iV ■ • • r> (Ais) 

|#>:=JV- 1/2 Xy <fci |$) (A16) 



becomes the first excited state with the gap Ai = a' . 
In contrast to the open XXZ model discussed in ap- 
pendix A 1, the gap for the thermodynamic limit is hence 
finite. Consequently, the coherence decay is always ex- 
ponential in this model as described in the main text. 



FromiJefflV'i) 
the eigenequation 



a+a'\„i,j\ a-a' n./.j-l 



Wi 



a + a a — a 



wr }+\^i )) follows 



cos(fc) |^ fc ). (A17) 



The minimum eigenvalue a' is obtained for k = 0. 

Using exact diagonalization (ED), we have asserted 
that the first excited state is similar to \<f>\) [Eq. (A13)] 
for a' /a ~ 1 (small J z ) with a gap Ai ~ a+ 4 a . For suffi- 
ciently small a' /a (sufficiently large J z ), |V>i) [Eq. (A16)] 



Appendix B: One and two spins with L = S z 

For completeness and in order to have a consistency 
check for the the perturbative analysis presented above in 
appendix A 2, let us shortly discuss the decoherence for a 
single spin-1/2 in a transverse magnetic field. The Liou- 
ville master equation for the Hamiltonian H = —2h x S x 
and the Lindblad operator L — S z reads 



dtp = i[2h x S x , p] + ~f(S z pS z - p/4). 



(Bl) 



Parameterizing the density, matrix as p = [ b ^ ic y2^ ] , 
we have 



d t p = ih x 



-2ic 1 -2a 

2a - 1 2zc 



6 + ic 
6-ic 



giving the three equations 



<9ta = 2h x c, dtb = 6, 9tc = h x (l — 2a) c. 

So 6 decays as b(t) = o(0)e^*/ 2 . The first and third 
equations imply dfc = — 4h x c — \dtc with the solution 
c{t) = c+e iw+t + c-e luj - f + c.c. and 



-=w^-a 



(B2) 



For 7 < 8 1 h x |, the coherence hence decays in the long- 
time limit as ~ e~ 7 */ 4 . For large 7 however (7 > 8|/i x |), 
the square root in Eq. (B2) becomes imaginary and, 
thus, alters the decoherence rate from 7/4 to 7/4 — 

a/(7/ 4 ) 2 ~ Ah l = %KH + £>(K/l 3 )- This can be com- 
pared with the perturbative upper bound (a + o/)/4 
[Eq. (A14)] for the gap of the transverse-field Ising model 
as derived in appendix A 2. For J z = 0, this bound has 
indeed the same value found here, (a + a')/ A = 8h x /j. 
So the two computations are consistent. 

In the main text, we also mention the decoherence for 
two spins with zero total magnetization (Sf + S*| = 0) 
evolving according to the Lindblad master equation with 
the XXZ interaction H' = ^(S+S^ + S^S}) + J Z S(S% 
and the dissipative term V'p' = ^'(Sfp'Sf + S%p'S% ~ 
///2). This model maps to the single spin in a transverse 
field (Bl) if we associate the state \\\) with the spin-up 
state, and |4-t) with the spin-down state. The resulting 
single-spin dynamics is governed by the Hamiltonian H = 



JxyS x — </z/4 and the dissipative term T>p = 2^'{S Z pS z — 
p/4:). So the equation of motion is just Eq. (Bl) with 
transverse field h x = —J xy /2 and bath coupling 7 = 27'. 
Hence, for 7' < 2|J X1/ |, the decoherence rate is 7'/2. As 
is shown in the main part of the paper, the decoherence 
changes from this exponential behavior to a power-law 
decay in the thermodynamic limit of the corresponding 
many-particle system. 



Appendix C: Convergence of the tDMRG results 

The numerical investigations are based on the 
time-dependent density matrix renormalization group 
(tDMRG) [26-28]. We use a Trotter-Suzuki decompo- 
sition of the propagator exp(£i) [34] for which the error 
is of fifth order in the time-step At = 0.125. As described 
in the main text, an essential part of the algorithm is to 
truncate small-weight components in the singular value 
decomposition of the density matrix p(t) . This is neces- 
sary in order to bound the computation cost. The corre- 
sponding truncation threshold e determines the accuracy 
and the computation cost of the simulation. One has 
to ensure convergence of the numerical results with re- 
spect to e and with respect to the system size N, as we 
are interested in the physics of the thermodynamic limit 
N -)• 00. 

Figure 5 shows the analysis for the dissipative XXZ 
model (Al). The results are sufficiently converged for 
a truncation threshold of e — 10~ 12 and the system size 
N = 96. The simulations presented in the main text were 
done with these parameters. Figure 6 shows the analysis 
for the dissipative transverse-field Ising model (A6). A 
truncation threshold of e = 1CP 10 and the system size 
N = 64 are sufficient and were used for the simulations 
presented in the main text. 
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